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Abstract 

We investigate the temporal evolution of impulsively generated pertur- 
bations in a potential coronal arcade with an embedded loop. As the initial 
configuration we consider a coronal loop, represented by a density enhance- 
ment, which is unbounded in the ignorable direction of the arcade. The 
linearized time-dependent magnetohydrodynamic equations have been nu- 
merically solved in field-aligned coordinates and the time evolution of the 
initial perturbations has been studied in the zero-/3 approximation. For prop- 
agation constrained to the plane of the arcade, the considered initial pertur- 
bations do not excite trapped modes of the system. This weakness of the 
model is overcome by the inclusion of wave propagation in the ignorable 
direction. The inclusion of perpendicular propagation produces two main 
results. First, damping by wave leakage is less efficient because the loop is 
able to act as a wave trap of vertical oscillations. Second, the consideration 
of an inhomogeneous corona enables the resonant damping of vertical os- 
cillations and the energy transfer from the interior of the loop to the external 
coronal medium. 

Subject headings: Sun: atmosphere - Sun: corona - Sun: oscillations - 
Sun: surface mangetism 



1. Introduction 



In the last years a phenomenon that has attracted the attention of solar physicists is 
the discovery of transverse oscillations of coronal loops observed in EUV wavelengths 
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(171 A) with TRACE in 1998. The oscillatory amplitude and period are of the order of 
a few Mm and a few min, respectively, and one of the most interesting features of these 
oscillations is that their amplitude decreases quickly with time, typically in a few peri- 
ods. Some examples of this oscillatory phenomenon, n amely flare-excited tran sversal 



o scillations, were reported by lAsc hwande n et al 
recent ly by 
(l201ll) . see 



Wans & Solanki (2004); Hor i et al 



Aschwanden et al. 



( 1999f): lNakariakov et all (I 199%. more 
J2005h. and lAsch wanden & Schriiver 



(Eq02D; iSchrijver et al.1 (120021) and lAschwanden! (120091) 
for an extensive overview and analysis of transversal flare-excited coronal loop oscilla- 
tions. 

The first theoretical studies of the modes of oscillation of coronal flux tubes mod- 



eled as strai ght magnetic cylinde rs were done by IWentzell (Il979() . ISpruid (| 198 lh . Edwin & Roberts 

(|1983|) and Roberts et al.l (| 19841) . Later, flare-generated transverse oscil lations were in- 

terpreted as fast kink eigenmode oscillations of straight c ylindrical tubes dNakariakov & Ofman 



200ll:lRuderman & Roberts 2002; 



nal ma gnetic structures, 



Goossens et al. 



Goossens et al 



2002). In the context of curved coro- 



(1985): IPoedts et al.l(ll985l) : IPoedts & Goossens 



(1988 ) investigated the c ontin uous spectrum of id eal magnetohydrodynamics (MHD). 
Oliver et al] (I1993L Il996|) and iTerradas et al] (|1999|) derived the spectrum of modes in 
potential and non-potential arcades. More complex configurati ons, such as sheared m ag- 
netic arcades in the zero plasma-/? limit, have been studied by lArregui et al] (|2004al ib|). 
Other authors have studied eigenmodes in curved configurations with density enhance- 
ments that represent coronal loop s (e.g., Van Doorsselaere et al. 20041; Terrada s et al 



2006b; Verwichte et al.ll2006al]b 



Diaz et al]l2006l : lvan Doorsselaere et al.ll2009h . The 



fact that the corona is a highly inhomogeneous and structured medium complicates the 
theoretical description of MHD waves and it is believed that this could be the underlying 
cause of the observed wave damping. Several mecha nisms of wave damp ing have been 



sorption ( 


Hollwesr & Yandl988l:lGoossensll991 


Goossens et al. 2002: 


Ruderman & Roberts 


2002; Van Doorsselaere et al] 


2004, and references therein) and more recently wave 


leakage by tunneling ( 


Bradv & ArbeJl2005: Bradv et al.ll2006: Verwichte et al.ll2006a b ct 



Diaz et al.l l2006). 



Although normal modes should be seen as the building blocks to interpret coronal 
loop oscillations they do not represent the whole picture, but their study provides a basis 
for understanding the dynamics of the system. To have a more accu rate description, the 
time-de pendent problem needs to be analyzed. Using this method, ICadez & Ballestei 
(|1995aUbl) studied analytically the propag ation of fast waves i n a two-dimensional coro- 
nal arcade with uniform Alfven speed. lOliver et al.l (|1998|) s tudied the effect of im - 
pulsively generated fast waves in the same coronal structure. Del Zanna et al] (|2005f) 
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studied the properties of Alfven waves in an arcad e configuration , includ ing a transition 
region between the photosphere and the corona. iTerradas et al.l (I2008b() used a poten- 
tial arcade embedded in a low environment to study the properties of linear waves. 
Other studies have analyzed the effect of the loop structure on the properties of fast and 
slow waves in two-dimensional curved c onfigurations (see , e.g., Brady & Arberl l2005 l : 



Terradasl (2005 



Murawski et alJl2005t iBradv et al.ll2006l : ISelwa et alJl2006L l2007h : see 
for a review. 

The aim of this work is to examine two physical mechanisms involved in the 
fast attenuation of the o bser ved vertical coronal loop oscillations. Previous studies by 
Verwichte et al. J2006affl : ISelwa et al.l Q2006, 2007) have considered the damping of 
vertical oscillations by wave leakage in a coronal loop embedded in a curved config- 
uration. When propagation is constrained to the plane of the structure, these authors 
find that wave leakage produces a strong damping and that the ratio of the damping 
time to the period is, in some cases, much shorter than the observed values. On the other 
hand, damping by resonant absorption, caused by the inhomogeneity of the medium, has 
mainly been studied in single m agnetic slabs (ITerradas et all 2005 : Arregui et al. 2007 ) 



Terradas et al 



2()()6at7for 



and in single magnetic cylinders iRuderman & Robert J2OO2I : 
example. Nevertheless, more complex equilibrium models have also been c onsidered 
dVan Doorsselaere et all2004l : ITerradas et al]l2006b[ l2008ah . iRial et all (bold) have re- 
cently considered the coupling of fast and Alfven modes in a potential coronal arcade 
with three-dimensional propagation of perturbations. This study shows that because of 
the inclusion of perpendicular propagation fast wave energy can easily be converted into 
Alfven wave energy at given magnetic s urfaces by means of resonant coupling. Our 
analysis aims at extending the model by [Rial et al.l (120101) by including a density en- 
hancement in a curved magnetic configuration in order to study how three-dimensional 
propagation affects the efficiency of the damping of vertical loop oscillations by wave 
leakage and how the inhomogeneity of the corona can produce coupling of modes and 
energy transfer. 

The paper is arranged as follows. In §|2]we describe the equilibrium configuration 
as well as the approximations made in this work. In §[3] we present the linear ideal MHD 
wave equations with three-dimensional propagation of perturbations. In §@]we describe 
the numerical setup together with the initial and boundary conditions. Our results are 
shown in § [51 where the linear wave propagation properties of coupled fast and Alfven 
waves in a two-dimensional coronal loop, allowing three-dimensional propagation, are 
described. Finally, in §|6]the conclusions are drawn. 
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2. Equilibrium configuration 



The equilibri um magnetic field is a potential arcade contained in the xz-plane (see 
Oliver et al.ll 1993L for more details). In Cartesian coordinates the flux function is 



A(.v. ;:) = B A t! cos ( J exp f--^- |. ( 1 ) 



and the magnetic field components are given by 



B x (x, z) = B cos |-^J exp J , 

B z {x, z) = -B sin |^-J exp ( -^-j . (2) 

In these expressions A B is the magnetic scale height, which is related to the lateral extent 
of the arcade, 2L, by A s = 2L/n, and B Q is the magnetic field strength at the base of the 
corona (z = 0). The overall shape of the arcade is shown in Figure [1^. 

In this paper gravity is neglected and the B = approximation is used. Under these 
assumptions the equilibrium density, p, can be chosen arbitrarily. We consider a loop 
with uniform density, p , embedded in a corona whose density, p e , is also uniform and 
smaller than that of the loop by a factor 10, i.e., p e = po/10; see Figured^. The vertical 
density profile at the arcade center is shown in Figure[2l 

The combination of the magnetic field of Equation © with this sharp density profile 
leads to the following Alfven speed distribution 

f v AO exp(-^-), inside the loop, 
v A (x,z) = \ } 7< . (3) 

I v Ae expy--^), otherwise, 

where v A o = B / yjpo/uo and v Ae = B / ^p e po are the Alfven speed inside and outside the 
loop at the base of the corona (z = 0). This formula gives v A at any point in the jcz-plane. 
Notice that the Alfven speed varies both along and across magnetic field lines in our 
curved configuration. 



Oliver et al.1 (|1996|) and lArregui et al.1 (|2004a|) showed that an appropriate descrip 



tion of normal modes in a curved structure, such as the one considered here, can be 
obtained by solving the MHD equations in flux coordinates, which are determined by 
the previous election of the equilibrium field in Equation ©. Appropriate flux coordi- 
nates are given by the following expressions 
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Fig. 1 . — (a) Magnetostatic configuration of the potential coronal arcade described by 
Equations (QQ) and ©, where the solid curves represent magnetic field lines, given by 
A(x, z) =constant. These curves in the xz-plane become arcade surfaces in three dimen- 
sions. In this model z measures the vertical distance from the base of the corona (placed 
at z = 0). The black region represents a coronal loop embedded in the arcade. Two 
orthogonal unit vectors defining the normal and the parallel directions, e n and e\\, are 
also shown at a particular point, (b) Same as (a) in the^-plane. Curved magnetic field 
lines in the xz-plane become straight lines in this flux coordinate system and therefore a 
curved loop becomes a straight slab. The four boundaries are color-coded to show their 
correspondence with lines and points in panel (a) 



«K*,z) = cos^— jexp^-— j, < <A < 1, (4) 

, . Sm (As) eX P( A B ) - 1 

x{x > z) = (1-^1/2 ' -1<*<1, (5) 

where iff andx are the coordinates across and along the equilibrium magnetic field lines. 
The normalization of the length of all field lines to the same value is achieved by the 
factor (1 - t/r 2 ) 1/2 in the denominator of Equation ©. One of the advantages of using 
this coordinate system is that it enables to include the whole coronal arcade in the finite 
domain ^ e [0, 1],^ e [-1, 1]. In this domain magnetic field lines are straight and each 
of them is represented by a different value of iff. In addition, the magnetic field strength 
depends on both if/ and^. The shape of the potential arcade and the coronal loop in these 
field related coordinates is shown in Figure Q}). 
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Fig. 2. — Vertical variation along the z-axis of the density (solid line) and Alfven speed 
(dotted line) for the coronal loop configuration of Figured] 



3. Magnetohydrodynamic equations and linear waves 

In order to study small-amplitude oscillations in our potential arcade with an em- 
bedded loop the previous equilibrium is perturbed. For linear and adiabatic MHD per- 
turbations in the zero-/? approximation the relevant equations are 



p_l = _(VxB 1 )xB, 
at no 



(6) 



= V x ( Vl x B), 



(7) 



where p and B are the equilibrium density and magnetic field and the subscript 1 is used 
to represent the perturbed velocity, Vi, and magnetic field, B^ 

When these equations are particularized to our two-dimensional equilibrium it turns 
out that it is advantageous to use field-related components instead of Cartesian compo- 
nents in order to characterize the directions of interest related to the polarization of each 
wave type. The unit vectors in the directions normal, perpendicular, and parallel to the 
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equilibrium magnetic field are given by 



VA 
|VA|' 



e i = e 



3" 
|B|' 



(8) 



where A is the flux function, given by Equation and |B| = B Q exp (-z/A B ) is the mag- 
netic field strength. In a low-/? plasma and in the absence of perpendicular propagation, 
these three directions are associated with the velocity perturbation of the three types of 
waves that can be excited, namely v\ n for fast waves, v lx for Alfven waves, and viu for 
slow waves. 

As the equilibrium is invariant in the y-direction, we can Fourier analyze all per- 
turbed quantities in the y-direction by making them proportional to exp (ik y y), where k y 
is the perpendicular wave number. In this way, three-dimensional propagation is allowed 
and each Fourier component can be studied separately. By making the appropriate def- 
initions, the field-related components of the MHD wave equations can be cast in the 
following manner 



dt 

dVly 

~dT 

dB ln 
dt 

dB iy 

dt 

dt 



|B| 

Hop 
IBI 



(B a - if/ a d„ -x a d\\)B n + (iff s d» +xA\ ~ B s) B \n 
(i// s d n +x.A\) B iy + k y B i 



= \B\(i// s d n +Xsd\\+B s )v ln , 
= \R\(i// s d n +x s d\\)v ly , 



-IBI 



kyViy + {ift a d n +x a d\\ + B a )v u 



(9) 

(10) 

(11) 
(12) 
(13) 



The notation used to designate the different components of the fields and coordi- 
nates is related to the directions defined by the equilibrium field, so that the field compo- 
nents with subscript n, y, and || refer to the directions normal, perpendicular, and parallel 
to the equilibrium field lines. To unify the notation, the partial derivatives related to the 
flux coordinates if/ andx have been also rewritten with the subscript n and || respectively, 
i.e. 
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8 = e ■ V 

8y =iky = e y -V, (14) 

d|| = h ■ V. 

Furthermore, the quantities ij/ s , i// a ,Xs,Xa, B s , and B a , that represent the derivatives 
of the flux coordinates and the magnetic field s trength along (subs cript Y) and across 



(subcript 'a') magnetic field lines, are defined in [Oliver et al.1 (|1996|) . 



fa = £(VA • iff, = J(B • V)^, 

*« = g(VA.V)*, Ar, = s(B-v)Ar, (15) 

5, = ^(VA • V)5, 5,. = £(B • V)fl. 

Notice that such a general derivation and implementation of the governing equa- 
tions enables us to apply the scheme to any other configuration for which similar flux 
coordinates can be defined. 

Equations (l9l)-(fT3l) constitute a set of coupled partial differential equations with 
non-constant coefficients that describe the propagation of fast and Alfven waves. As the 
plasma 8 is zero, slow waves are excluded from the analysis and viy = 0. When k y = 0, 
Equations (l9l)-(fT3l) constitute two independent sets of equations. The two equations for 
v\ y and B\ y are associated to Alfven wave propagation. On the other hand, the three 
equations for the remaining variables, v\ n , B\ n , B\\\, describe fast wave propagation. 



4. Numerical method and initial and boundary conditions 



The obtained set of differential equations are too complicated to have analytical or 
simple numer ical solutions. Fo r this reason we have solved them with a numerical code 
described in lBona et al.l (120091) . using flux coordinates. This particular choice of coordi- 
nates is crucial for the proper computation of the solutions, and has been overlooked in 
previous numerical studies. 

We try to mimic the observed vertical oscillations of coronal loops in the corona 
when a sudden release of energy occurs and perturbs the structure. Since we are con- 
cerned with the excitation of fast waves, our initial perturbation is such that only the nor- 
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mal velocity component is disturbed, whereas all the other variables (vi y , B\ n , B\ y , B\\\) 
are initially set to zero. To initially excite the system we have chosen a two-dimensional 
profile in V\„ with, respectively, a cosine function and a Gaussian profile along and across 
field lines, 



This expression represents an initial disturbance that perturbs a range of field lines cen- 
tered about the field line ip = if/ . Here vo is the initial amplitude of the disturbance, a 
is related to the range of field lines initially affected by the perturbation, and k\\ is the 
wave number of the disturbance along the magnetic field. This symmetric initial profile 
does not represent a general disturbance in the solar corona, but it is selected in order to 
mainly excite the fundamental fast mode with one maximum along field lines. This can 
be achieved by an adequate choice of k\\. A more general disturbance would excite an 
ensemble of modes at the same time, making the interpretation of the obtained numerical 
results more difficult. 

The implementation of the appropriate boundary conditions in the numerical code 
is an important issue. The reflection of waves at the bottom boundary, due to the large 
inertia of the photospheric plasma, is accomplished by imposing line-tying boundary 
conditions at z = 0. In all other boundaries, flow-through conditions are imposed so 
that perturbations are free to leave the system. As the equations are solved in flux co- 
ordinates, one then needs to make the proper translation between the system boundaries 
from Cartesian to field related coordinates; see Figure [0 Although the whole arcade can 
be reproduced in flux coordinates, considering the complete range < ifr < 1 causes nu- 
merical issues, so we restrict ourselves to the range [0.034, 0.9] in the i/r-direction. This 
implies that extremely low and extremely high field lines are discarded in the numerical 
simulations. 



The results presented in this section have been obtained with the numerical solution 
of Equations (l9l)-(fT3l) after an initial perturbation given by Equation (fT6l) with ifr = 
0.41 and a = 0.05 is launched. Different values of k y have been considered (k y L = 
0, 5, 16, 60). The two-dimensional variation of v\„ and v\ y for some of these simulations 
is presented as animations associated to Figures [3] and [51 Time in these animations and 
in subsequent plots is given in units of ta = L/vao- 




(16) 



5. Numerical results 



- 10- 




-1.00 -0.50 0.00 0.50 1.00 

Fig. 3. — Snapshot of the two-dimensional distribution of the normal velocity compo- 
nent for ky = 0. Some magnetic field lines (black curves) and the edge of the coronal 
loop (white lines) have been represented. (An animation of this figure is available in the 
online journal.) 



5.1. Wave propagation in the plane of the arcade — wave leakage 

We first consider wave propagation in the plane of the arcade (k y = 0). Such as can 
be appreciated in the animation accompanying Figure[3l the initial perturbation produces 
traveling disturbances across the magnetic surfaces that propagate away from the density 
enhancement. After a short time (J/ta ~ 5) the loop contains very little energy, and so 
these disturbances can be interpreted as a combination of leaky modes that the loop 
structure is unable to confine. This interpretation is also supported by the time evolution 
of the normal velocity component at the center of the loop (Figure Sk), that displays a 
strong damping such that, for tlx a ^ 5, the velocity has an almost null amplitude inside 
the loop. 

A convenient way to quantify the wave energy leakage and the damping is to use 
the total energy density computed both in the full domain and in the interior of the loop. 
We use the following formula to calculate the energy density at a given position and 
time, 
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Fig. 4. — Results for k y = 0. (a) Temporal evolution of the normal velocity component, 
v ln , at the loop center (x = 0, z/L = 0.56). (b) Normalized energy density in the whole 
domain (solid thick line), inside the loop (dotted line), and outside the loop (dashed line) 
as a function of time. The two-dimensional temporal evolution of v\ n is presented in the 
animation associated to Figure [3] 



The total energy density in a spatial domain, D, of the xz-plane can be computed from 
the spatial integration of 6E(r, t), 



Jd 

In our plots this quantity is normalized by dividing it by the initial energy density over 
the whole numerical domain, E(0). 

Figure|4}) shows the total energy density integrated over the full domain, the interior 
of the loop, and the external region. We can see that the loop is unable to trap wave 
energy and that, after a little more than 5ta, it has already transferred all its energy 
to the outside medium. After this time, the wave energy outside the loop equals the 
energy in the full domain. Because of this wave leakage, the total energy of the system 
continuously decreases and when the last leaky waves carrying a non-negligible amount 
of energy reach the domain boundaries (i.e. at tlx a ~ 10), it is practically zero; see also 
the animation of Figure [3] 



6E(t, t) = \ p [v\ n + v\ y ) + - (B 2 ln + B\ y + fl? ) . 

Z /in 



(17) 




(18) 
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These re sults differ from wh at is obtained in a line-tied straight slab model of a 
coronal loop (Terradas et al.ll2005|) . in which the kink mode of the loop remains oscillat- 
ing after the initial leaky phase, since the density enhancement is able t o trap part of the 



energy of the initial perturbation . Such as was already po inte d out by iMurawski et al. 



(120051) : ISelwa et all (I2005L 120061) : IVerwichte et al.l (l2006aUbUcl) : ISelwa et all (120071) . the 
magnetic field line curvature produces damping by wave leakage and the energy de- 
posited initially in the loop is simply radiated as a combination of leaky modes, and thus 
no trapped solutions are found. 



5.2. Wave propagation in three-dimensions — wave trapping 



iTerradas" 



et al.l (|2006b) showed that in a toroidal loop struct ure wave leakage is 



much less effective, that in slab geometry. Following the same idea. lArregui et al.l(|2007|) 
used a simple three-dimensional straight slab model to demonstrate that in the case of a 
slab, introducing oblique propagation increase s the spatial confinement of the kink mode 



around the loop. IVan Doorsselaere et al.l(|2009) have also shown that a cylindrical tube is 



a more efficient wave guide than the slab model, and less energy is allowed to leak. When 
a potential arcade is used as the equilibrium structure, the oblique propagation of the slab 
corresponds to perpendicular propagation, which means the introduction of the k y wave 
number. If the previous results for the kink mode with oblique propagation in a coronal 
slab are also valid in the present loop configuration, we expect that an initial perturbation 
will excite both leaky waves (that will leave the system in a similar manner as found in 
§ 15.11) and trapped modes (that will retain energy in the system). A comment about the 
term "trapped" mode used in this paper is in order. In our plasma configuration the local 
cut-off frequency depends on position because of the Alfven speed inhomogeneity and 
the inclusion of perpendicular wave propagation. For this reason, an eigenfunction that 
is evanescent in a certain region can change its character to propagating in another part 
of the system, and so wave energy leakage can arise. For this reason we warn that modes 
termed "trapped" here may actually loose energy by this mechanism. If leakage is very 
small the modes have, to all practical purposes, a trapped behaviour. 

As in § 15.11 we have considered an impulsive excitation of the normal velocity 
component, but now k y L = 5. The temporal evolution of the normal and perpendicular 
velocity components (see the animations accompanying the top panel of Figure [5]) illus- 
trate that there are important differences compared to the case k y = (Figure [3]). First 
of all, some of the energy contained in the initial perturbation goes to the perpendicu- 
lar variables because k y ± implies that the perturbed perpendicular components (y\ y 
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k y L = 5 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

x/L x/L 



-1.00 -0.50 0.00 0.50 1.00 -1.20 -0.60 0.00 0.60 1.20 

k„L = 60 




-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 

x/L x/L 



-1.00 -0.50 0.00 0.50 1 00 -1.00 -0.50 0.00 0.50 1.00 

Fig. 5. — Snapshots of the two-dimensional distribution of the normal and perpendicular 
velocity components for k y L = 5 and k y L = 60. Some magnetic field lines (black curves) 
and the edge of the coronal loop (white lines) have been represented. (An animation of 
this figure is available in the online journal.) 



and Biy) are no longer independent of the normal and parallel ones (vi„, B\ n , and B\\\). 
In addition, we have just discussed that for k y = most energy has left the system at 
tlx a ~ 10. For k y L = 5, however, it is apparent that some energy remains in the loop 
and part of the arcade above it for much longer times. This is interpreted as a signature 
of the excitation of one or more trapped modes of the system. We next substantiate this 
claim with a careful inspection of v\„. 

We consider the temporal variation of the normal velocity component near the loop 
center (Figured). It displays two distinct behaviors: for t/r A < 5 very short periods 
are prominent. They correspond to leaky modes, that propagate quickly away from 
the loop. After this leaky phase all that remains is a gentle, damped oscillation. This 
damping will be later discussed in detail, but let us mention that it is milder than that 
of the k y = case (Figure Hk). The power spectrum of the V\ n signal (Figure (6fc) shows 
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two clear peaks that, according to the interpretation in the previous paragraph, simply 
reflect the power contained in two trapped modes excited by the initial perturbation. An 
evidence of this would be the presence of power peaks at nearly the same frequency 
in neighboring points. To investigate this possibility, we consider v\ n as a function of 
t for x = and different values of z, compute their power spectra and stack them in a 
contour plot (Figure[6fc). The result is the presence of significant power at two horizontal 
contours whose frequencies coincide with those found in Figure [6}), each of them with 
its characteristic range of heights. From Figure |6fc we see that the fundamental mode 
frequency is coL/vaq ~ 0.95 and that this mode covers the range 0.5 < z/L < 1.15 for 
x = 0. The first harmonic has a>L/v A0 ~ 1.28 and covers the range 0.3 < z/L < 0.95 for 
x = 0. 

Hence, we conclude that two trapped modes of the system have been excited by the 
initial perturbation and that this is the reason why the normal velocity component is less 
strongly attenuated for k y = 5 than for k y = 0. Both modes are spatially distributed over 
the coronal loop and a large region above it, that is, perpendicular wave propagation 
increases the wave confinement in the present equilibrium configuration, although for 
k y L = 5 wave energy around the loop is also found. 




Fig. 6. — Results for k y L = 5. (a) Temporal evolution of v in near the loop center (x = 0, 
z/L = 0.52). (b) Normalized power spectrum of the signal of panel (a), (c) Power 
spectrum of V\ n at x = as a function of height. The three solid lines are, from bottom 
to top, the frequency of the fundamental Alfven mode and its first two harmonics. The 
two dashed lines mark the limits of the coronal loop. 
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5.3. Mode coupling and resonant energy transfer 

Additional information about the temporal evolution of the system comes from 
the perpendicular velocity component. In Figure [7_k this perturbed variable is plotted 
at two different positions at the arcade center (x = 0). In both cases v ly displays an 
oscillatory behavior with a monotonically increasing amplitude and it is evident that the 
oscillatory period is different at the two positions. This is confirmed by the power spectra 
of these two signals, each displaying a single peak at one of the frequencies of the normal 
velocity components found in Figures [6b and c. Next, power spectra of vi y as a function 
of t along the z-axis are computed and stacked together (Figure |7fc) and it becomes 
clear that most of the power of this velocity component is concentrated at two particular 
heights. To explain this result we recall that in § 15 .21 we mentioned that after the initial 
perturbation, trapped modes of the system were established by a transfer of energy from 
v\ n to the other perturbed variables, including v\ y . But a second, more efficient, process 
takes place here: the resonant transfer of energy from the trapped modes to the Alfven 
modes whose frequencies match those of the former. Figure [6]: shows that the frequency 
of the two prominent trapped modes found in § l5.2l (q)L/v4n ~ 0.95 and u)L/v A o ~ 1.28) 
intersect the lowest one of the three solid lines included in this plot. These solid lines are 
the frequencies of th e Alfven continua corresponding to the fundamental mode and its 
first two harmonics ([Oliver et al. 1993 ). Hence, we conclude that the frequencies of the 



two trapped modes match that of the fundamental Alfven mode at two particular heights, 
namely z/L ~ 1.12 and z/L ~ 0.96, so that these are the positions where Alfven energy 
appears, in agreement with Figure |7fc. 

It should be noted that the trapped mode frequencies also match the frequency of 
Alfven harm onics, but this happens for heights at which the trapped mode amplitude 



is negligible. lArregui et al.1 (I2004a|) showed that for k y ± resonant coupling between 
fast and Alfven modes can only happen for modes with the same parity along B. As 
a consequence, resonant absorption with Alfven modes other than the fundamental one 
does not take place. 

So far our description of the resonant absorption process is incomplete because it 
relies on the information along the z-axis. The animation of the top panel of Figure [5] 
shows that the largest amplitudes of vi_ v for t/r A > 10 are in two ranges of magnetic 
surfaces centered about the field lines whose apexes are at z/L ~ 0.96 and z/L ~ 1.12 
(i.e. at the resonant positions in Figure [7];). A relevant issue is that the Alfven frequency 
varies in each of these two ranges of magnetic surfaces and so the resulting Alfven 
oscillations soon become out of phase. This phase mixing creates small spatial scales 
that lead to the numerical dissipation of energy whose importance will next become 
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Fig. 7. — Results for k y L = 5. (a) Temporal evolution of v iy at the points x = 0, 
z/L = 1.12 (solid line) and x = 0, z/L = 0.96 (dotted line), (b) Normalized power 
spectrum of the signals of panel (a), (c) Power spectrum of V\ y at x = as a function of 
height. The three solid lines are, from bottom to top, the frequency of the fundamental 
Alfven mode and its first two harmonics. The two dashed lines mark the limits of the 
coronal loop. 

clear. 

Equation (fT8l gives a measure of the energy density in a domain of our system from 
the spatial integration of Equation (flTT ). In these definitions, 8E(v,t) and E(t) contain 
the contribution from both the Alfvenic perturbed variables (i.e. v ly and Bi y ) and the 
perturbed variables characteristic of fast modes for k y = (i.e. v\ n , B hu and B\\\). We 
denote by E^ifvenif) and Ef ast (t) the spatial integral of these two contributions over the 
whole system, 

E A ifven(t) = [pv\ y + ±B\ y ) dxdz, 

(19) 

E fast (t)= S\[pv[ n + ±(B\ n + B^}dxdz. 

Since these two integrals are carried out over the full arcade, it is clear that E M f ven {t) + 
E/ast(t) is equal to E(t) of Equation (fT8l) when D is the complete numerical domain. 

Figure [8k shows the energy density in the system, inside the loop and outside it. 
The first of these three quantities presents three phases: for t/r A < 5 the total energy 
in the whole system decreases rapidly because of the emission of leaky waves. Then, 
for 5 < tlx A < 25 the rate of energy decrease is slower, and for t/r A > 25 it becomes 
stronger again. To understand the last two phases one must bear in mind the energy loss 
by the numerical dissipation of Alfven waves. During the phase 5 < t/r A < 25 there 
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is a substantial amount of energy in the loop (dashed line in Figured) that is slowly 
transferred to the surrounding medium (dotted line in Figure [8^). We know that this 
process is caused by resonant absorption, by which the energy that could not leak in 
the initial phase is transferred to the two resonant positions. After t/r A ~ 25 almost 
all the energy in the trapped modes has been given to Alfven modes, which implies 
that these modes cannot longer "feed" from their previous energy "reservoir". As a 
result, the numerical dissipation of phase-mixed Alfven waves proceeds at a faster pace. 
Figure [8J5 provides support to this interpretation: at t/r A ~ 10 almost all the energy has 
been converted from the "fast"-like perturbations to the Alfvenic ones, but Ef ast does 
not become negligible until t/r A ~ 30. This means that some conversion to E A if ven still 
takes place for 10 < t/r A < 30, and so the dissipation of Alfven waves is slower than in 
the later phase. 




Fig. 8. — Results for k y L = 5. (a) Normalized energy density integrated over the whole 
system (solid line), integrated over the coronal loop (dotted line), and integrated over 
the region outside the loop (dashed line), (b) Normalized energy density integrated over 
the whole system contained in all perturbed variables (solid line; same as solid line of 
panel (a)), in the "fast" perturbed variables, Ef ast (i) (dotted line), and in the Alfvenic 
perturbed variables, E A[ f Ven (t) (dashed line). 
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5.4. Wave confinement and resonant absorption for large k y 

To evaluate the influence of the perpendicular wave number on the confinement 
of wave energy by trapped modes of the loop and the loss of this energy by resonant 
absorption to Alfven continuum modes we consider k y L = 60. We first inspect the 
two-dimensional temporal evolution of the perturbed velocity components (animations 
associated to the bottom panel of Figure [5]). There are obvious differences with the 
numerical simulation with k y L = 5: first of all, there is a lack of wave leakage at the 
beginning of the simulation, shortly after the initial impulse is released. Second, large 
velocity amplitudes are only found inside the loop structure, both for the v ln and the v ly 
components. And third, since the excited trapped modes are not spread outside the loop, 
there is no resonant absorption at large heights. In summary, these animations seem 
to indicate a much stronger confinement of the wave energy by the loop. Next, let us 
confirm this preliminary result. 

Now, we consider the temporal variation of the normal and perpendicular velocity 
components along the z-axis and compute the power spectra of these two signals, which 
are then stacked to produce a contour plot of the power as a function of height along 
the z-axis and frequency. Figure [9] reveals that the power of v\ n an d V\ y is concentrated 
inside the loop and that these signals contain two frequencies around coL/v A o = 0.6. 
A comparison of these graphs with Figures [6t and |7J; gives a clear confirmation of the 
stronger wave confinement achieved by the loop for k y L = 60, or in other words, the 
stronger spatial confinement of the trapped mode of the system. Moreover, since the 
trapped mode does not extend to large heights, resonant absorption is irrelevant in the 
present numerical simulation. For resonant absorption to be possible, the trapped mode 
should have power at a height z/L ~ 1.5. 

We finally pay attention to the energetics. Figure [10b shows that during the whole 
simulation there is almost no energy outside the loop, in agreement with our previous 
findings about the wave guiding properties of the loop for a large perpendicular wave 
number. Regarding the energy transfer from the "fast" to the Alfven perturbed variables, 
Figure [TOb indicates that while some energy is transferred from the first to the second 
one, on the long term the "fast" components retain a larger proportion of energy. The 
total energy of the system presents a decay in time that, unlike to that described in § 15.31 
is constant during the whole simulation. There are two reasons for this difference: first, 
the loop is such a good wave guide for transverse oscillations when k y L = 60 that very 
little wave energy is carried out of the system by leaky waves. Then, the initial leaky 
phase of Figure [8] is not found in Figure \\0\ Second, for k y L = 5 resonant absorption 
brings energy to Alfvenic motions at large heights, where the spatial mesh is coarser. 
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Fig. 9. — Results for k y L = 60. Power spectrum of (a) V\ n and (b) v ly at x = as a 
function of height. The three solid lines are, from bottom to top, the frequency of the 
fundamental Alfven mode and its first two harmonics. The two dashed lines mark the 
limits of the coronal loop. 



This enhances numerical dissipation of Alfven waves, which results in a strong wave 
damping. For k y L = 60, however, all wave energy is retained at small heights, where the 
spatial mesh is denser. Hence, although numerical dissipation cannot be removed from 
our numerical experiments, it is greatly reduced. 



5.5. Realistic three-dimensional coronal loop 



Some authors (see for example iHollweg & Yangj|l988|) have extended the results 
of the slab models to the cylindrical geometry by using the equivalent of the azimuthal 
wave number in Cartesian coordinates. The formula that relates both wave numbers is 



k y = -, (20) 
r 

where r is the loop radius and m is the azimuthal wave number. Since we are interested 
in vertical oscillations, caused by the kink mode, we take m = 1. Regarding the loop 
radius, it is obvious from Figure [Ik that it changes with height, being largest at the 
apex and smallest at the feet. Since the physical conditions at the loop top are the most 
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Fig. 10. — Results for k y L = 60. (a) Normalized energy density integrated over the 
whole system (solid line), integrated over the coronal loop (dashed line), and integrated 
over the region outside the loop (dotted line), (b) Normalized energy density integrated 
over the whole system contained in all perturbed variables (solid line; same as solid line 
of panel (a)), in the "fast" perturbed variables, Ef ast (t) (dotted line), and in the Alfvenic 
perturbed variables, E A if ven {t) (dashed line). 



determinant for the kink mode properties, r is taken as the loop half width at the apex. 
Then we obtain k y L =16. 

To study the properties of vertical oscillations for this value of k y we first consider 
the temporal evolution of V\ n at a point inside the loop (see Figure [TTa). The periodic 
behaviour that was found for k y L = 5 is also present in this case, although the attenua- 
tion rate is smaller now (compare with Figured). This suggests a better confinement of 
the vertical oscillations for k y L =16 and to quantify this effect we compute the power 
spectra of V\ n along the z-axis and plot them together as a contour plot, which is pre- 
sented in Figure [Tib. It is clear that two trapped modes, with frequencies ojL/vao = 0.74 
and coL/vao = 0.88, are excited by the initial perturbation and that they are spatially 
spread along a wide range of heights. Although most of the energy is contained in the 
fundamental trapped mode and, in particular, inside the loop, the two trapped mode fre- 
quencies match the fundamental continuum Alfven frequency around z/L = 1.1-1.3. 
This opens the possibility of resonant absorption playing a role in the damping of the 
trapped modes. To investigate this effect, we plot the perpendicular velocity component 
at two points along the z-axis (Figure [TTb). 
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The first of these two signals is taken inside the loop, at its top boundary, while the 
second signal is gathered at one of the resonant positions. Their amplitude aside, these 
two signals are different in that the first one has an appreciable amplitude just after the 
initial leaky phase (i.e. for t/r A > 5), whereas the second one is initially negligible and 
only after t/r A ~ 10 starts to grow, at a rate larger than that of its counterpart inside 
the loop. To understand these two behaviors we plot the power spectrum of vi y as a 
function of height and frequency (cf. Figure [TTb). which confirms the presence of large 
power in the transverse velocity component both inside the loop and in the resonant 
position. The power in the range z/L = 1.1 - 1.3 is concentrated around the Alfven 
continuum frequency and we interpret this as the signature of resonant absorption. On 
the other hand, inside the loop (i.e. in the range z/L = 0.5 - 0.64 there is power at the 
Alfven continuum frequency but also for coL/v A q ~ 0.74, that is, at the frequency of the 
fundamental mode detected in the normal velocity component. We thus conclude that 
the perpendicular velocity component inside the loop exists both as part of the trapped 
mode excited by the initial perturbation and by the resonant transfer of energy from this 
trapped mode to Alfven continuum modes. 



6. Conclusions 

In this paper we have studied the temporal evolution of fast and Alfven waves in 
a curved coronal loop embedded in a magnetic potential arcade, in order to asses the 
relevance of three-dimensional propagation of perturbations on the damping of vertical 
loop oscillations by wave leakage and resonant absorption. 

When perpendicular propagation is not included (i.e. when waves are constrained 
to propagate in the plane of the arcade), a transverse impulsive perturbation produces a 
combination of leaky modes and the loop is unable to trap energy in the form of vertical 
kink oscillations. The energy deposited initially in the loop is emitted rapidly to the 
external medium. This result confirms previous findings in the sense that in a curved 
coronal loop slab model damping by wave leakage is an efficient mechanism for the 
attenuation of vertical loop oscillations. 

When three dimensional propagation of waves is considered, two new effects are 
found. First, damping by wave leakage is less efficient and the loop is able to trap part of 
the energy deposited by the initial disturbance in transverse oscillations. The amount of 
energy trapped by the structure increases for increasing values of the perpendicular wave 
number. Second, the inclusion of the perpendicular wave number in a inhomogeneous 
corona produces the resonant coupling between fast and Alfven modes at those posi- 
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Fig. 1 1 . — Results for k y L = 16. (a) Temporal evolution of the v ln velocity component at 
x = 0, z/L = 0.61. (b) Temporal evolution of the vi v velocity component at x = 0, z/L = 
0.61 (solid line) and at x = 0, z/L = 1.2 (dotted line). Power of the normal velocity 
component for x = as a function of z and the dimensionless frequency, (d) Power of 
the perpendicular velocity component for x = as a function of z and the dimensionless 
frequency. In panels (c) and (d) the solid l ines are the theoretical frequencies of the 
Alfven continua given by lOliver et al.1 (|1993|) . 



tions where the trapped mode frequency matches that of Alfven waves. This resonant 
coupling produces the transfer of energy from the fast wave components to Alfvenic 
oscillations. In our model the loop boundary is sharp and so there is no smooth density 
transition that allows resonant absorption to happen in this position. Hence, this energy 
transfer does not occur at the loop boundary, but at locations in the external medium, in 
particular above the coronal loop. 

Let us stress that slab models of curved coronal loops in the absence of perpendic- 
ular propagation give rise, in general, to damping times by wave leakage shorter than 
those observed. On the contrary, as we have demonstrated in this work, perpendicular 
propagation is the clue to obtain damping times compatible with observations. 
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